Matter X waves 
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We predict that an ultra-cold Bose gas in an optical lattice can give rise to a new form of 
condensation, namely matter X waves. These are non-spreading 3D wave-packets which reflect 
the symmetry of the Laplacian with a negative effective mass along the lattice direction, and are 
allowed to exist in the absence of any trapping potential even in the limit of non-interacting atoms. 
This result has also strong implications for optical propagation in periodic structures. 
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Matter waves are a natural manifestation of large scale 
coherence of an ensemble of atoms populating a funda- 
mental quantum state. The observation of Bose-Einstein 
condensates (BECs) in dilute ultra-cold alkalis [l| has 
initiated the exploration of many intriguing properties of 
matter waves, whose mascroscopic behavior can be suc- 
cessfully described via mean-field approach in terms of 
a single complex wave-function with well defined phase 
across the atom cloud 2]. Large scale coherence effects 
are usually observed by means of 3D magnetic or optical 
confining potentials (also ID cigar-like or 2D disk-shaped 
BEC are possible 3]) in which BECs are described by 
their ground-state wave-function However, trapping can 
also occur in free space (i.e. without a trap) through the 
mutual compensation of the leading-order (two-body) in- 
teraction potential and kinetic energy, leading to bright 
(dark) solitons for negative (positive) scattering lengths. 
This phenomenon has been observed only in ID (i.e., soli- 
ton waveguides, a trap confining the BEC transversally 
to the soliton) 4]. In 2D and 3D, free-space localization 
cannot occur due to collapse instability of solitons, and 
even in a trap collapse usually prevents stable formation 
of BEC 0, 0| , needing stabilizing mechanisms . 

A lot of attention was also devoted to periodic po- 
tentials due to optical lattices 0, where the behavior of 
atom mimic those of electrons in crystals or photons in 
periodic media 9|, and exhibit effects which stem from 
genuine atom coherence In ID (elongated) lattices 
bright (gap) solitons can form also in the presence of re- 
pulsive interactions 111. In this letter we predict that 
a novel trapping phenomenon occurs when the full 3D 
dynamics is retained in a ID lattice. Specifically, un- 
der the conditions for which the Bloch state associated 
with the lattice has a negative effective mass, the natural 
state of the BEC is a localized matter X wave character- 
ized by a peculiar bi-conical shape. The atoms are or- 
ganized in this way in the absence of any trap, solely as 
the result of the strong anisotropy between the ID mod- 
ulation and the free-motion in the 2D transverse plane. 
Furthermore, due to axial symmetry, the atoms can ex- 
perience a collective motion with given velocity along the 
lattice, resulting into wave-packets traveling undistorted. 
Remarkably, these matter X waves are not only observ- 



able individually, but allow also to describe physically 
realizable BECs as their superposition, each component 
maintaining a constant (in time) number of atoms. 

The very first step of our derivation is similar to that 
of matter gap solitons. Starting from the (mean-field) 
Gross-Pitaevskii (CP) equation [3| with an optical stand- 
ing wave potential and no additional trap (we set r] = 



ihdtijj 



-r/V^V7 + 4rsin2(fcz/2)V' + a|V'|V = 0, (1) 



we assume axial symmetry around z and decompose the 
wave-function ip — ip(r,z,t) (r^ = -{- y^) into its for- 
ward and backward components as 



(2) 



which, for slowly varying in z, and dropping rapidly 
rotating terms, allows us to reduce the CP equation (Q) 
to the coupled equations 



(3) 



where £± = ihdt±irikdz+ri'V\, and = df.+r^^dr- In 
the linear limit a = 0, the plane-wave [exp(iKz — iEt/fi)] 
linear dispersion relation associated with Eqs. Q has 
two branches E — E±(k) — ±.Ty^l + (we set p = 
K.krj/T), exhibiting an energy gap of width 2r. The 
coupling between il)f and -0;, causes the structure to be 
strongly dispersive near band-edge and the linear dynam- 
ics of atoms to be governed by strong Bragg reflection. 
Nevertheless, in the ID limit (ry = 0), where Eqs. © were 
obtained previously ,lj|, the nonlinearity (both attrac- 
tive a < and repulsive a > 0) induces self-transparency 
mediated by a two-parameter family of moving bright gap 
solitons, so-called because they exist in the gap seen in 
the soliton moving frame In the attractive case, one 
might think that the nonlinearity can balance also the ki- 
netic transverse term Vj_ leading to bell-shaped 3D atom 
wave-packets |0. We show in the following that, con- 
trary to this expectation, close to the lower band edge 
E ~ E^ , the atomic wave-function takes a completely 
different form. To this end we apply a standard enve- 
lope function (or effective mass [lij |) approximation |l4l |. 
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searching for spinor solutions tp — [tpf of the form 

^ = e0(er, ez, et)~^ _ exp(?Kz - itE^/fi) + O(e^) (4) 

where e is a small expansion parameter^ is slowly modu- 
lating the Bloch state with amplitude V - — [V'/- i'b-V 
(eigenvector of Eqs. (O with a = corresponding to 
the eigenvalue E-). At the leading order we find that 
obeys the following asymptotic equation 

ifidtcj) + iE'_d^(l) + E'Ldlct) + ryVi^ - xl'/'P^ = (5) 

where x = and ^ ^-^^ E'L ^ ^ account 

for dispersion. For sake of simplicity we deal henceforth 
with the strict band-edge case k — Q, which can be pre- 
pared by acting on the wave-number and on the potential 
parameter In this case Eq. Q reads explicitly as 

where —m^ = —mT/{r]k'^) is the negative effective mass 
associated with the lattice, in turn determining the hy- 
perbolic character of (GP) Eq. The scaling transfor- 
mation z,r,t,(j> — > ZQZ,rQr,tot, Co4> with Tq = ZQUie/m, 
to = 2meZQ/h, Cq = 2?i/(3|a|io), being a length scale, 
allows us to use dimensionless variables. In the attractive 
case (a < 0) , Eq. JHJ) admits matter X waves solutions of 
the kind (j) = ip{r, z) exp{i^t), where ip{r, z) is indeed an 
X-shaped invariant envelope 0|. Notice that X waves, 
well known in optics 0, acoustics or microwaves 
HOl, as non-spreading (in space and time) solutions of 
the linear Helmotz wave equation j2lj , have been only re- 
cently discovered for Schodinger-type models of the form 
of Eq. lO, by analyzing the so-called paraxial (linear 
[2^ or nonlinear 0|) propagation in dispersive media. 
However, we emphasize that, while an optical or acous- 
tic field retains a directly observable spatio-temporal X- 
shape 0,01, in the case of a Bose gas, the local density 
of atoms IV'P has an X-shaped spatial envelope mod- 
ulating a term cos^(fcz/2) (periodic with lattice period), 
due to the form of Eq. (O . 

To assess further the regimes of observability of matter 
X waves, we discuss two crucial issues. First, by extend- 
ing the analysis of Ref. 0| , we are able to show that X 
wave solutions of Eq. JB)) exist also in the (far more com- 
mon) case of repulsive nonlinearity (a > 0). As an exam- 
ple, we display in Fig. ^Jthe atom density corresponding 
to the stationary solution with eigenvalue n = 0, which 
shows a dense core accompanied by bi-conically shaped 
regions of lower density. As shown in the insets, the sig- 
nature of the X-shape is a single peak on-axis (r — 0) 
and a double peak off-axis (r ^ 0). The envelope |0p 
exhibits also slow oscillations and modulate the fast si- 
nusoidal variation of the density IV'P- 

Second, our aim is to show that non-spreading atom 
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FIG. 1: (Color online) Surfaces of constant (high, interme- 
diate, and low in blue, red, and bronze, respectively) enve- 
lope atom density \cj>\'^ of a matter envelope X wave solution 
(/I = 0) of Eq. @. The insets show on-axis {r = 0) and off- 
axis (r — 7) longitudinal profiles of the overall , black) 
and envelope (|</>P, red) densities. We fix zo = W/k to have 
dimensionless units. 

X-shaped BECs can be formed also in an ideal noninter- 
acting gas (a — 0), which is realizable by exploiting Fesh- 
bach resonance |^ . This is in contrast with other known 
settings where BEG needs either a confining potential 
(as in pioneering experiments ^]^), or the nonlinearities 
to balance kinetic spreading. Ultimately, this stems from 
the fact that X waves have a finite linear limit, unlike soli- 
tons of standard (elliptic) GP equation whose amplitude 
vanishes as a ^ 0. To demonstrate, however, that linear 
matter X waves are observable, we need to address the 
fundamental issue of their norm finiteness. To this end 
let us solve Eq. (jSJ in dimensionless form with a = 

{tdt + vi - af) cf, = o, (7) 

by seeking for envelopes moving with velocity v 

(j){r, z, t) = ip{r, C) exp(-TO2:/2 -I- itv'^/4), (8) 
where = z — vt, and ip turns out to obey the equation 
d^if + r-^drip - dl^p = 0. (9) 
The general solution of Eq. (jSJ 

/>oo 

^{rX) = / /(a)Jo(ar)e'"^da, (10) 
Jo 

represents a class of envelope X waves specified by their 
spectrum f{a), which generalize to (Schrodinger-t ype ) 
Eq. {Tj) the X wave solutions of Helmotz equation [21|. 
An exponentially decaying (with arbitrary inverse width 
A) spectrum /^(a) = exp(— aA), yields the simplest (or 
fundamental) X wave ipx = p -I- (A-<)2] while 
/i""* = a"exp(— aA) (n = 1,2, ...) defines the derivative 
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where F{a,v) = ct~^ J2q^'}i''^)fq~i'^)- introducing 
new variables kt, kz such that a = kt and v = 2{kt — kz), 
and setting U{kt,kz) = 2F{kt,2kt - 2kz), Eq. JHI) can 
be cast in the form 



»s = 



X waves = d^'ifx/dA" . In general from Eq. ((TUI) . we The importance of Eq. (fTK|l stems from the fact that 

obtain the full envelope X wave solution of Eq. as (f>s describes a wide class of physical atom beams. To 

^ show this, we start from the integral representation of 

0(r, z,t)= I /(a) Jo(ar)e*(""^)"+*('^""'')*da, (11) 4>q in term of its spectrum CSJ, which yields 

^ poo POO 2 

where \^\^ clearly travels undistorted along z. However, ^ I I F{a,v)Jo{ar)e'^"-^-^'^-''^'adadv, (17) 

we have to face the pitfall that these X waves do not 
represent physical objects since their norm diverges. This 
follows from the transverse scalar product < (j>\(j> >± (z, t) 
[isj which yields for any pair of solutions (/), of Eq. 
with different spectra /, / (but equal velocity v) 

/>oo 

<(j)\(j) >j_^2Tr f{a)f*(a)a-^da. (12) 

^0 

Since < (f)\(f) >± docs not depend on z and t, the 3D 
norm < (f)\(f) > of any X wave (f) diverges, thus requiring 
an (unphysical) infinite number of atoms. Remarkably, 
however, finite norm beams can be generally constructed 
by introducing new orthogonal X- waves |24| . Inspired by 
Eq. H12|l . we can exploit the orthogonality of associated 
Laguerre polynomials Lg^'' (x) {q — 0, 1,2, ...) with respect 
to the function x exp(— x), to introduce a numerable class 
of (transversally) orthogonal X waves cj)'^ {r, z,t\v) defined 
by the following spectra (and parametrically by their ve- 
locity v) 



oo ^0 



Uikukz)Jo{ktr)e"'''e'^''"~''"^'ktdktdkz, 



7rv/2(9+l) 



iW(2Aa)e- 



-Aq 



(13) 



which represents the generic axisymmetric solution of Eq. 
Q, expressed in 3D momentum space (kt,kz), kt being 
the momentum transverse to the lattice direction z. Here 
U{kt,kz) is the Fourier-Bessel (or plane- wave) spectrum 
of the initial atom distribution (f)o(r,z) = 4>{r,z,t = 0), 
and U{kt,kz) = 2k-^J2qCqi^kt - 2kz)f^ikt) stands 
for its expansion in terms generalized Laguerre poly- 
nomials [any square-integrable f{x) in x € [0, oo) can 
be expanded in terms of l'^\x)]. This argument can 
be reversed by stating that, given the initial (t = 0) 
distribution of atoms U{kt,kz) in momentum space, if 
F{a,v) = U{a,a — v/2)/2 is square integrable with re- 
spect to a (with a € [0, oo)), then the atom wave-packet 
admits the representation IjlSfl . The expansion coefficient 
can be easily calculated as 



The X waves cj)-^ satisfy the orthogonality relation < 
l<^^ >j_= Spq/An, with Spq the Kronecker symbol. Im- 
portantly, the 3D scalar product shows that such waves 
are orthogonal also with respect to the velocity, i.e. any 
pair of waves with velocity u and v satisfies the relation 

< (l)^{r,z,t\v)\(j)j{r,z,t\u) >^dpqd{v~u). (14) 

From Eq. (|14|l it is natural to consider the solution (j) = 
of Eq. given by the superposition of orthogonal 
X waves {4>q} travelling with different velocities v as 

/oo 
Cq{v)(p^{r,z,t\v)dv. (15) 

q -o° 

From the orthogonality relation (|14|l . we find that the 
total number of atoms is 



Cq{v) 



V2A 



= [ U (a,a--) e-2"^iW (2aA)ada, 
1 Jo ^ 2/ 



(16) 



where N'q = \Cq{v)\'^dv represents the atom number 
of the q-th X wave component 4''^{v) of the wave-packet. 
Therefore we obtain the remarkable result that, while the 
superposition (p-^ generally describes atom wave-packets 
which evolve in time, such evolution preserves the distri- 
bution of atom number among the X wave components. 



Clearly, most of physically relevant wave-packets belong 
to the class (f>Y, 26]. For example a spectrally narrow 
gaussian beam can be described by few X waves (details 
will be given elsewhere). 

Once the existence of finite norm linear X waves is 
established together with their potential to describe gen- 
eral BECs, the most intriguing question remains whether 
we can expect matter X-shaped atom distributions to 
be observable. Such waves correspond to a single fixed 
value q = q in Eq. (|15|l and an ideal velocity distribution 
Cq{v) — S{v — v). However, more generally, we can con- 
sider an atomic envelope beam (j)a constituted by several 
replicas of the single X wave (f>^ travelling with different 
velocities. Although, strictly speaking, such beam is not 
stationary, it can approximate such a state with an ar- 
bitrary degree of accuracy. In other words, it is possible 
to construct solutions that preserve their shape, for an 
arbitrary long time. Indeed if Cq{v) is a narrow function, 
e.g. peaked around v = 0, Eq. H15I) . using also Eq. (jS)), 
yields the following atom envelope 0q 



Mr,z,t) = / ^fir,z-vt)Cqiv)e~'^''+''^'dv 



ipUr, z)c{z,t) 



(18) 
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FIG. 2: Atom density of a pure X wave {(pti^y with g = 
and A = 1, as seen on the {x,z) plane {y — 0). The insets 
show (same units as in Fig. 1) the overall (thin line) and 
envelope (bold line) density on-axis and off-axis, respectively. 



where c{z,t) = J^^Cq{v) exp{—i^z + i^t)dv is a so- 
lution of the dispersive wave equation idtc — die = 0. 
Eq. (fH^ represents an X wave modulated by a dispers- 
ing wave, which gives a sort of adiabatic dynamics of the 
finite norm X wave. Indeed the atom beam has an in- 
variant spatial shape fixed by (f^, which we display in 
Fig. [2 as an example for q — 0. It spreads on a charac- 
teristic time which is longer, the narrower is the velocity 
distribution function Cq{v). Finally, while in the linear 
regime we expect that such atom states should be some- 
how prepared, we envisage that atom collisions (nonlin- 
ear regime) can strongly favour the formation of X waves 
from more conventional ball-shaped atom clouds (e.g., 
obtained by a harmonic 3D trap which is then switched 
off) through instability mechanisms 0, H^, an issue 
which will be deepened elsewhere. 

In conclusion we have shown that a periodic poten- 
tial supports moving or still localized states of the GP 
model with envelope X-shape (strictly or nearly) pre- 
served upon motion. A matter X waves entails local- 
ization both in momentum and configuration space and 
is a clear signature of a Bose condensed gas, so much as 
the anisotropy in the distribution function detected 
in early experiments. However, unlike any other form of 
BEG including solitons, matter X waves can be observed 
in free-space and in the non-interacting regime, where 
they constitute a natural basis of expansion to describe 
the coherent properties of atom wave-packets. These re- 
sults have strong implications also in optics, where the 
model (O holds for normally dispersive bulk media [j^ . 
or in nonlinear optics of stratified media where Eqs. (jSJ 
describe propagation in a ID bulk grating in the pres- 
ence of diffraction, or propagation in a fully 3D photonic 



crystal (along proper directions). 
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